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ABSTRACT 



The pulsating DA white dwarfs are the coolest degenerate stars that undergo self- 
driven oscillations. Understanding their interior structure will help to understand the 
previous evolution of the star. To this end, we report the analysis of more than 200 h of 
time-resolved CCD photometry of the pulsating DA white dwarf star EC 14012—1446 
acquired during four observing epochs in three different years, including a coordinated 
three-site campaign. A total of 19 independent frequencies in the star's light variations 
together with 148 combination signals up to fifth order could be detected. We are 
unable to obtain the period spacing of the normal modes and therefore a mass estimate 
of the star, but we infer a fairly short rotation period of 0.61 ± 0.03 d, assuming the 
rotationally split modes are i = 1. The pulsation modes of the star undergo amplitude 
and frequency variations, in the sense that modes with higher radial overtone show 
more pronounced variability and that amplitude changes are always accompanied by 
frequency variations. Most of the second-order combination frequencies detected have 
amplitudes that are a function of their parent mode amplitudes, but we found a few 
cases of possible resonantly excited modes. We point out the complications in the 
analysis and interpretation of data sets of pulsating white dwarfs that are affected by 
combination frequencies of the form /a + /s — fc intruding into the frequency range 
of the independent modes. 

Key words: stars: variables: other ~ stars: white dwarfs - stars: oscillations - stars: 
individual: EC 14012—1446 - techniques: photometric 



1 INTRODUCTION 



White dwarf stars are the most common end point of stel- 
lar evolution. Almost all stars with masses below SM© end 
their lives in the white dwarf stage, and most of them do so 
after being the central star of a Planetary Nebula. Once the 
central star has exhausted its nuclear fuel, it slowly cools 
and dims at nearly constant radius. As there is a definite 
low-temperature cutofi' to the white dwarf luminosity func- 
tion, it can be used to determine the age of the galactic disk 
(e.g. Winget et al. 1987, Richer et al. 2000) or stellar clusters 
(e.g. Hansen et al. 2007), although some theoretical uncer- 



tainties remain (see Fontaine, Brassard & Bergeron 2001 for 
a recent review). 

While white dwarf stars cool, they cross three pulsa- 
tional instability strips, i.e. regions where they exhibit self- 
driven nonradial oscillations. The hottest instability region 
contains pulsating PG 1159 stars (the GW Vir stars), the 
intermediate one consists of pulsating DB white dwarfs (the 
V777 Her stars), whereas pulsating DA stars (the ZZ Ceti 
stars) constitute the coolest white dwarf instability strip. 
If all white dwarfs showed self-excited pulsations as they 
pass through their instability domain, this would mean that 
the interior structure of the pulsators is representative of 
all white dwarfs. The ZZ Ceti instability strip is considered 
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to be pure by most authors (e.g. see Bergeron et al. 2004, 
Mukadam et al. 2004, Castanheira et al. 2007). 

In any case, if a large number of oscillations can be 
detected and resolved in a given star, one can apply aster- 
oseismic methods to it. In other words, the oscillations can 
be used to determine the interior structure of such stars. 
This is important for astrophysics in general, as the history 
of evolution of stars is engraved in the interiors of its end 
product, the white dwarf star. Fortunately, white dwarfs are 
well suited for asteroseismology because virtually complete 
pulsational mode spectra over a certain frequency interval 
were found in a number of cases (e.g. Winget et al. 1991, 
1994), for instance for the prototype pulsating PG 1159 and 
DB white dwarf stars. 

The pulsating DA white dwarfs have remained less ac- 
cessible to this method. Their interior structure is somewhat 
more complicated due to their surface hydrogen layer and 
possible interior crystallization (e.g. Montgomery & Winget 
1999). Their pulsational mode spectra are in general less 
dense than those of the pulsating PG 1159 and DB stars, 
and even so, some of the observed signals are difficult to 
be reconciled with theoretical models (e.g. Pech, Vauclair & 
Dolez 2006, or Castanheira et al. 2004). 

Part of these difficulties may arise from combination fre- 
quencies, i.e. signals that may only reflect the Fourier decom- 
position of the non-sinusoidal light curve shapes and that do 
not correspond to normal mode oscillations, intruding into 
the frequency domain of the normal modes. Another poten- 
tial problem is the occurrence of amplitude and frequency 
variations on short time scales (e.g. Handler et al. 2003) , re- 
sulting in spurious peaks that complicate the identification 
of normal modes. On the other hand, the temporal varia- 
tions of the pulsational spectra of some ZZ Ceti stars may 
make modes that are not permanently observable visible at 
some periods of observation. This fact was first recognized 
and taken advantage of by Kleinman et al. (1998). 

Still, considerable progress in the understanding of ZZ 
Ceti star pulsation has been made recently. For instance, 
Kepler et al. (2005) pointed out how the C/O ratio in the 
core of these stars can be measured from their evolutionary 
period changes that refiect their cooling rate. Furthermore, 
based on the original idea by Brickhill (1992), Montgomery 
(2005) derived a method to constrain pulsational mode iden- 
tifications from the light curve shapes of pulsating white 
dwarfs, and to recover the thermal response time scale of 
the convection zone, which depends on effective tempera- 
ture. 

Consequently, it is desirable to obtain extensive obser- 
vations of multiperiodic ZZ Ceti stars, and to decipher their 
pulsational mode spectra. A good target for asteroseismic 
studies must be carefully chosen: in general, hot pulsating 
DA white dwarfs have few modes, low amplitudes and short 
periods, whereas the cooler pulsators have longer periods, 
higher amplitudes and more modes, but temporal variations 
in their pulsation spectra. The best choice may then lie in 
between: when observing pulsators in the centre of the ZZ 
Ceti instability strip, short and long periods, high ampli- 
tudes, as well as stable mode spectra may be present. 



2 OBSERVATIONS AND REDUCTIONS 

The ZZ Ceti star EC 14012-1446 was discovered to be 
a multiperiodic high-amplitude pulsator by Stobie et al. 
(1995), who found five independent and five combination 
frequencies in their data. As this object seemed well suited 
for asteroseismic purposes, we have obtained time-resolved 
CCD photometry of the star over four epochs distributed 
over three years. Three of the data sets were acquired at a 
single site, but one data set was obtained during a coordi- 
nated campaign involving three observatories well spread in 
geographical longitude. 

The first measurements originated from the 0.75-m 
telescope at the South African Astronomical Observatory 
(SAAO) and were taken in April 2004. A high-speed CCD 
photometer (O'Donoghue 1995) was used in frame-transfer 
mode with integration times of 10 seconds. No filter was 
employed. 

Two months later, we observed EC 14012-1446 with 
the 1.9-m telescope at the SAAO with the same photometer. 
Due to the smaller field of view, frame-transfer mode was 
not employed, but the integration time of 10 seconds was 
kept, resulting in one frame per 12 seconds, and again no 
photometric filter was inserted. 

In the year 2005, we organized a multi-site campaign 
for the star, involving three observatories in South Africa, 
the United States and Taiwan. In this way, close-to uninter- 
rupted measurements could be obtained during the central 
part of the campaign. The 1.0-m telescope at Lulin Obser- 
vatory started the observations, which were taken as unfil- 
tered 12-second integrations, resulting in one frame per 17 
seconds. The 1.0-m telescope at SAAO came on line four 
nights afterwards, this time using a standard CCD with 21- 
second exposures and no filter, yielding one data frame every 
30 seconds. Later during the same night, the 0.9-m telescope 
of the Southeastern Association for Research in Astronomy 
(SARA) located in Arizona joined the campaign, acquiring 
unfiltered measurements with 12-second exposures for one 
data point per 16 seconds. 

Finally, in 2007, the 0.9-m telescope at the Cerro Tololo 
Interamerican Observatory (CTIO) was used to measure 
EC 14012—1446 again. These data were acquired through 
a S8612 red-cutoff filter. The integration time was 10 sec- 
onds, leading to a data rate of one frame per 30 seconds. 
The star was observed every second night of this run. An 
overview of all the observing runs is given in Table 1. 

All measurements from the years 2004 and 2005 were 
reduced in the same way. First, standard IRAF0 routines 
were used to correct the images for overscan, bias level (if 
needed), dark counts (if needed) and flat field. Photome- 
try was carried out using the MOMF (Multi-Object Multi- 
Frame, Kjeldsen & Frandsen 1992) package. MOMF applies 
combined Point-Spread Function/ Aperture photometry rel- 
ative to an optimal sample of comparison stars, ensuring 
highest-quality differential light curves of the target. 

The CTIO data from 2007 were also reduced with IRAF 
in the above-mentioned way, but photometry was carried out 
with a series of IRAF scripts employing aperture photom- 

^ IRAF, the Image Reduction and Analysis Facility, is written 
and supported by the IRAF programming group at the National 
Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
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Table 1. Overview of the observations. AT is the total time of 
monitoring, and Af is the frequency resolution of the data set. For 
the multisite campaign, the frequency resolution of the combined 
data is quoted and marked with an asterisk. 



1 - 



Month/ Year 


Telescope 


# nights 


AT 


Af 








(hr) 


(mHz) 


April 2004 


SAAO 0.75-m 


4 


36.4 


2.2 


June 2004 


SAAO 1.9-m 


4 


18.6 


3.7 


May 2005 


Lulin 1.0-m 


7 


26.8 


0.75* 


May 2005 


SAAO 1.0-m 


8 


50.9 


0.75* 


May 2005 


SARA 0.9-m 


5 


30.5 


0.75* 


April 2007 


CTIO 0.9-m 


7 


40.1 


0.87 


Total 




35 


203.3 





etry optimized for high-speed CCD data (Kanaan, Kepler 
& Winget 2002). The final output from these scripts again 
is a differential light curve of the target star. Tests showed 
that both of the CCD photometry programs employed by 
us yield results of comparable quality. 

All differential light curves were visually inspected and 
data obtained during poor photometric conditions were re- 
moved. As EC 14012—1446 is substantially bluer than the 
nearby stars used as comparisons, some differential colour 
extinction present in the light curves was carefully removed 
with the Bouguer method (fitting a straight line to a plot 
of magnitude vs. air mass). No low- frequency filtering was 
performed to avoid aflfecting possible intrinsic long-period 
signals. Finally, as no periodicities shorter than 140 seconds 
were present in the light curves, some data were merged to 
speed up the computations to follow: the April 2004 data 
were co-added to 20-second bins, the June 2004 data to 24- 
second bins, the Lulin measurements from 2005 to 34-second 
bins and the SARA data from 2005 to 32-sccond bins. This 
procedure also ensured that the 2005 campaign data all have 
similar time sampling and therefore similar weight in fre- 
quency analyses. 

Finally, all measurements were transformed to a com- 
mon time base. Wo chose Terrestrial Time (TT) as our refer- 
ence for measurements on the Earth's surface and applied a 
correction to account for the Earth's motion around the so- 
lar system's barycentre. As this barycentric correction varied 
by up to —1.7 s during single observing nights, we applied it 
point by point. Our final time base therefore is Barycentric 
Julian Ephemeris Date (BJED). The resulting time series 
was subjected to frequency analysis. 



3 FREQUENCY ANALYSIS 

The frequency analysis was mainly performed with the pro- 
gram Period98 (Sperl 1998). This package applies single- 
frequency power spectrum analysis and simultaneous multi- 
frequency sine-wave fitting. It also includes advanced op- 
tions such as the calculation of optimal light-curve fits for 
multi-periodic signals including harmonic and combination 
frequencies. Our analysis will require these features. 

We began by computing the Fourier spectral window 
functions of all data sets, and show the results in Fig. 1. 
The multisite campaign (May 2005) resulted in a reason- 
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Figure 1. The spectral window functions of our data sets from 
the four different epochs. 



ably clean window, whereas the other data sets suffer from 
alieising. The April 2004 window has strong sidelobes not 
only at the daily alias 11.6;ttHz away from the correct fre- 
quency, but also at 2.8, 8.8 and 14.3 /iHz, respectively. The 
June 2004 measurements have the poorest window function 
with low frequency resolution and strong diurnal aliases, and 
the April 2007 data suffer not only from daily aliasing, but 
also from spurious signals 5.2, 6.4, 16.9 and 18.0 /iHz away 
from the correct frequency. This has to be kept in mind when 
carrying out a frequency search. 

We proceeded by computing the Fourier amplitude 
spectra of the individual data sets, which can be found on 
the right-hand side of Fig. 2. Both the light curves on the 
left-hand side of Fig. 2 as well as the amplitude spectra 
make it immediately obvious that the dominant time scale 
of the light variations of EC 14012—1446 increased and that 
the amplitudes of individual pulsational signals also changed 
over time. 

Under these circumstances (individual data sets with 
time bases that are short compared to their separations, 
often complicated window functions, variable periods and 
amplitudes) it becomes clear that the frequency analysis will 
be difficult and that the data sets cannot be analysed jointly. 
We therefore used a different strategy: initiating the analysis 
with the apparently best data subset, and then using the 
experience gained for examining the other subsets. 
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Figure 2. Some example light curves (left-hand side) and the amplitude spectra of our data sets from the four different epochs. The 

dotted vertical lines denote the average distance between consecutive light maxima in the April 2004 data; note the temporal increase of 
the "mean" variability period directly visible in the light curves. 



3.1 The May 2005 multisite data 

We started the analysis by computing Fourier amplitude 
spectra of the light curves, and then determined the fre- 
quencies of the strongest signals. Those frequencies served 
as input parameters for a fit to the light curves that was 
optimized by nonlinear least squares sine- wave fitting imple- 
mented in Period98. In case of ambiguities due to aliasing, 
the frequency that yielded the lowest residuals between light 
curve and fit after optimization was chosen as input. As the 
next step, we prewhitened the fit from the light curve, com- 
puted the Fourier amplitude spectrum of the residuals and 
used it to search for further periodic signals. 

After a few prewhitening steps it became clear that am- 
plitude and/or frequency variations had occurred during the 
multisite campaign. This imposes severe complications on 
frequency analyses with Fourier algorithms, whose funda- 
mental assumption is that the amplitudes and frequencies of 
the individual signals present are constant within the data 
set. We therefore decided to use only the central part of the 
data set for the purpose of frequency search. This core paxt 
of the run contained 69% of the campaign data, spanned 



5.6 d, corresponding to a frequency resolution of 2.1 /xHz 
and had a duty cycle of 53%. 

Using only the central part of the run, we carried on 
with the prewhitening process as described above. Each sig- 
nal detected during this procedure was checked whether or 
not it was a combination frequency, i.e. whether it corre- 
sponded to a harmonic, sum, or difi^erencc frequency of sig- 
nals detected previously, within the frequency resolution of 
the data set (Table 1). If so, its frequency was fixed to the 
exact sum or difference of its paxent signals within PeriodOS. 
This procedure is valid here because the changes in the com- 
bination frequencies traced the changes in the parent modes 
within the accuracy of our data. 

In some cases, particularly for possible combinations of 
higher order, ambiguities arose because more than one set of 
parent frequencies could be matched to a given combination 
within the resolution of the data set. We then chose these 
parent frequencies whose product of amplitudes was highest, 
and required that the parent fi-equencies must have higher 
amplitudes than the combination signal they generate. 

Under these assumptions it is possible that the frequen- 
cies of some independent modes would be misidentified as 
combinations. However, we believe that disregarding a possi- 
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ble mode frequency is the safer choice compared to possibly 
accept some signals as real modes that are not. 

We continued the analysis until no significant peak was 
left in the residual amplitude spectrum, adopting the signal- 
to-noise criterion by Breger et al. (1993, 1999). In brief, an 
independent signal must exceed an amplitude signal-to-noise 
of 4 to be accepted, whereas S/N > 3.5 is sufficient for 
combination signals. We refer to the two papers mentioned 
before for detailed discussions. In the end, 15 independent 
frequencies and 90 combination signals up to fifth order were 
detected. 

3.2 The April 2004 data 

Because of the long nightly time series acquired during our 

first observing run on the star (more than nine hours per 
night on average), this single-site data set is also well suited 
for frequency analyses. With the same procedures as de- 
scribed before, we detected 12 independent and 44 com- 
bination frequencies in this data set, with no evidence for 
amplitude and/or frequency variations throughout the ob- 
servations. 

3.3 The June 2004 data 

Although acquired with the largest telescope we had avail- 
able for this study, this data set is not ideal for frequency 
analysis. The spectral window function is complicated and 
the frequency resolution is poor (see Fig. 1). Only the 
strongest signals could be detected without ambiguity, and 
we had to resort to adopting frequencies determined from 
the April 2004 data to avoid aliasing problems when try- 
ing to reveal signals with poorer signal-to-noise. Our fre- 
quency solution for this data sot containing 10 independent 
and 27 combination frequencies is probably incomplete. On 
the other hand, we believe that we have resisted the temp- 
tation to over-interpret this data set. 

3.4 The April 2007 data 

These measurements have the best formal frequency resolu- 
tion, but due to their sampling the spectral window is com- 
plicated. Still, we could determine a large number of the 
signals in the light curves independently, and only in a few 
cases wo had to adopt frequencies known from analyses of 
the other data sets. Amplitude and/or frequency variations 
likely occurred during these observations as well, but not on 
a scale severely hampering our analysis. Finally, we arrived 
at 11 independent and 40 combination frequencies for this 
data set. 

3.5 Putting it all together 

After the analyses of the individual data sets were com- 
pleted, the results were merged and the frequency solutions 
were compared. We list all 19 independent frequencies in 
our light curves in Table 2, all 52 first-order combination 
frequencies in Table 3, all 60 third-order combination fre- 
quencies in Table 4, and all 29 fourth-order combination 
frequencies together with the 7 fifth-order combinations in 
Table 5. 



4 DISCUSSION 

4.1 The independent frequencies 

As can be seen in Table 2, some, but not all of the inde- 
pendent frequencies are consistently present in all data sets. 
Some are singlets, and some show (regular) frequency split- 
tings. This requires individual discussions of all these signals 
or signal groups. Again, we start with the simplest cases and 
move on to the more difficult ones, choosing to proceed from 
the highest frequencies to the lowest. 

The signal at 2504 /iHz is well-behaved. It is present 
throughout all data sets, with its frequency and amplitude 
roughly staying constant, and it shows no splittings. It is 
the best candidate to search for possible evolutionary period 
changes, but our data are not sufficient to reach a definite 
conclusion in this respect. We note, and caution, that in 
the discovery data by Stobie et al. (1995) this signal had an 
amplitude 60% higher than in our measurements. 

The 1887 /LtHz variation is reasonably stable as well. 
However, in most data sets it had a companion near 
1881 /uHz, which is identified with a third-order combination 
frequency. The amplitude of the 1887 /iHz signal increased 
with time. This finding must be taken with caution because 
there is no information about the star's behaviour between 
the difi'erent data sets. This signal also showed a higher am- 
plitude in the Stobie et al. (1995) data than in any of our 
data sets. 

A frequency triplet centred at 1633 /iHz is present in all 
data sets. This triplet is always equally spaced within the er- 
rors, but the amount of the splitting varied: from 10.06 /iHz 
in April 2004 to 9.73 jjRz in June 2004, to 9.59 /^Hz in May 
2005 and to 10.44 ^Hz in April 2007. Given the amount by 
which this splitting changes and the observational errors, 
this variation is probably intrinsic to the star. 

An independent frequency at 1474 /iHz is also detected 
in each individual data set. In May 2005 it was a singlet, in 
June 2004 and April 2007 it was accompanied by another 
independent frequency lower by 10.8 and 9.1 /iHz, respec- 
tively, and in April 2004 it even was the centre of a triplet 
split by 9.98 /iHz. Within the errors, this splitting is consis- 
tent with the separation of the triplet around 1633 /xHz in 
the same period of observation (April 2004). 

A strong signal at 1385 /iHz is present in all data sets. In 
April and June 2004 it had a companion at a frequency 9.97 
and 10.60 /iHz lower, respectively. However, in April 2007 
the 1385 ytiHz signal was the centre of an equally spaced fre- 
quency triplet separated by 14.48 juHz. In May 2005, the 
situation is even more difficult: we found complicated struc- 
ture in this frequency region, with altogether six peaks and 
with the presence of two triplets centred on 1385 /iHz with 
both the « 10 and w 14 /iHz splittings. We note that the 
separation of some of the neighbouring signals is close to the 
resolution of our data set (1.5/AT, where T is the time base 
of the data, Loumos & Deeming 1978), but these signals re- 
main at the same frequencies when the full campaign data 
set is considered. 

At the two later observing epochs, a strong signal at 
1300 /iHz emerged. It was the strongest peak the amplitude 
spectrum in April 2007 as a singlet, and in May 2005 it 
had two lower-frequency companions separated by 3.3 and 
9.7 /iHz, respectively. There was no trace of it in our mea^ 
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Table 2. Independent frequencies in our light curves of EC 14012—1446. The formal errors in the amplitudes are ±0.4 mmag for April 
2004, ±0.5 mmag for June 2004, ±0.3 mmag for May 2005 and ±0.5 mmag for April 2007 







April 2004 






June 2004 






May 2005 






April 2007 




ID 


Free}. 


Ariipl. 


Frcq. 


Ampl. 


Frcq. 


Aiiipl. 


Free] . 


Ariijjl. 




(jillz) 




(//liz) 


(uiiiiagj 


(/,llz) 


(iiiiiiagj 


(//liz) 


(liiiiiag) 


/i 


821.26 ± 0.06 


7.1 


821.52 ± 0.13 


7.8 














/2 


1132.89 


± 0.14 


2.9 




















/3 














1289.28 


± 0.06 


5.3 








U 














1295.73 


± 0.04 


8.7 








h 














1299.00 


± 0.01 


39.8 


1300.62 


± 0.01 


63.6 


fa 














1371.89 


± 0.04 


9.1 


1370.89 


± 0.01 


20.0 


fr 


1375.53 


± 0.08 


5.5 


1375.02 


± 0.15 


6.9 


1375.60 


± 0.07 


4.7 








h 














1381.97 


± 0.03 


9.6 








h 


1385.50 


± 0.02 


24.4 


1385.62 


± 0.03 


35.1 


1384.84 


± 0.01 


40.6 


1385.32 


± 0.01 


44.1 


fio 














1394.91 


± 0.04 


7.3 








hi 














1398.29 


± 0.03 


10.9 


1399.84 


± 0.02 


13.6 


/12 


1464.17 


± 0.13 


3.2 


1464.17 


± 0.15 


6.8 








1463.95 


± 0.04 


6.4 


/l3 


1474.12 


± 0.04 


9.3 


1474.95 


± 0.11 


9.1 


1473.02 


± 0.06 


5.6 


1473.04 


± 0.04 


6.3 


/14 


1484.13 


± 0.16 


2.6 




















/16 


1623.28 


± 0.04 


11.2 


1623.51 


± 0.11 


9.0 


1623.86 


± 0.03 


10.5 


1623.20 


± 0.01 


18.6 


/16 


1633.36 


± 0.01 


48.1 


1633.60 


± 0.02 


48.3 


1633.71 


± 0.01 


33.5 


1633.69 


± 0.05 


4.7 


fl7 


1643.40 


± 0.03 


14.2 


1642.96 


± 0.07 


13.6 


1643.03 


± 0.05 


7.0 


1644.08 


± 0.06 


4.1 


fis 


1887.47 


± 0.05 


8.9 


1887.79 


± 0.11 


9.3 


1887.34 


± 0.03 


12.3 


1887.59 


± 0.02 


15.2 


/l9 


2504.86 


± 0.05 


8.7 


2504.98 


± 0.12 


8.1 


2504.97 


± 0.05 


6.8 


2504.65 


± 0.03 


8.6 



surements from 2004, and it was not detected by Stobie et 
al. (1995). 

The two lowest-frequency independent variations in our 
light curves were only present during (part of) 2004: a sin- 
glet at 1133 /iHz in April, and another singlet at 821 /xHz 
in both observing periods in 2004. We note that although 
numerically consistent in June 2004, an interpretation of /i 
as the subharmonic of fir generated a poor fit to the data 
and was therefore rejected. 

Stobie et al. (1995) reported an additional independent 
frequency of 1067 /iHz in their discovery data. However, it 
was only present in their first night of observation, and we 
cannot confirm its presence in any of our data sets. As we 
do not have the data by Stobie et al. (1995) at our disposal, 
we will not discuss this signal any further. 

4.1.1 Short-term amplitude and frequency variations 

As noted before, we had to concentrate our frequency anal- 
ysis of the May 2005 multisite photometry on the central 
part with the highest duty cycle because amplitude varia- 
tions had occurred during the observations. After prewhiten- 
ing the multifrequency solution listed in Tables 2 — 5 from 
these data, we still find residual mounds of amplitude left 
around most of the independent frequencies, most notably 
in the region with the most complicated structure around 
1385 AiHz. 

This raises the suspicion that at least some of the sig- 
nals in this frequency domain are artefacts due to amplitude 
and/or frequency variability. However, tests to check this hy- 
pothesis remained inconclusive. Consequently, it is not clear 
whether or not all the signals around 1385 /iHz (Table 2) 
indeed correspond to normal mode frequencies. 

We examined the occurrence of amplitude and fre- 
quency variations in the May 2005 data. We used the full 



data set for this purpose, and selected six independent sig- 
nals, either singlets or the strongest multiplet components. 
First, all but these signals were prewliitcned from the light 
curves, using Tables 2 — 5 as input and assuming constant 
amplitude and phase for all variations. We subdivided the 
residual light curves into sections long enough to avoid beat- 
ing within the six independent signals under consideration 
and computed the amplitudes and phases for each. These 
are represented in graphical form in Fig. 3. 

It is hard to say something quantitatively about the 
amplitude and frequency changes (which would require an 
even higher observational duty cycle) , but we can make some 
qualitative statements. 

• All independent signals do show amplitude and fre- 
quency modulation, with the possible exception of /19. 

• In general, signals at lower frequency (and therefore 
higher radial overtone) show more erratic behaviour. 

• Amplitude variations arc always accompanied by phase 
(frequency) variations. Whenever the amplitude is constant, 
the frequency is also constant. 

Figure 3 also confirms that the amplitudes and phases of the 
variations remained constant during the central part of the 
observations, validating the most important assumption of 
the frequency analysis in Sect. 3.1. 

Finally, we point out the interesting behaviour of signal 
/i3 (top right panels in Fig. 3) : it showed two phase "jumps" 
of ~ TT rad, at the beginning and the end of the observations. 
This might be an efi^ect of beating between closely spaced 
signals. The amplitude spectrum of the full data set indeed 
shows two peaks of similar amplitude spaced by 1.3 /iHz at 
this frequency. However, such a close frequency doublet is 
not present at a significant level in the measurements from 
2007 that would also have the time base to resolve it. 
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Table 3. Second-order combination frequencies in our light curves of EC 14012—1446 





Apri 


1 21)0 1 


luiic 


200 1 


Ala\- 


2005 


April 


2007 


ID 


Frcq. 


Ampl. 


Frcq. 


Ampl. 


Frcq. 


Ampl. 


Frcq. 


Ampl. 




(mHz) 


(mmag) 


(mHz) 


(mmag) 


(mHz) 


(mmag) 


(a*Hz) 


(mmag) 


/e - /s 














70.27 


3.9 


/g - /s 










85.83 


5.9 


84.69 


6.7 


/l2 - /s 














163.33 


3.1 


/i6 - /g 


247.86 


2.8 


247.98 


4.9 


248.87 


7.1 






/i5 - /e 














252.31 


3.9 


/l8 — fl6 


254.12 


3.6 














fi7 - h 






267.94 


5.9 










/i5 — /s 














322.58 


2.9 


/l6 - /o 










334.71 


6.0 






/l7 - /s 














343.45 


3.0 


/i8 - /g 










502.50 


1.8 


502.27 


2.5 


/l8 — /s 










588.34 


1.9 


586.97 


5.0 


/ig - /i5 














881.44 


3.5 


/ig — h 














1133.75 


5.2 


fi+fr 


2196.79 


3.3 


2196.55 


5.4 










h+h 


2206.76 


2.4 














h + /l3 


2295.39 


1.6 


2296.48 


2.7 










/l + /l6 


2454.62 


3.0 


2455.13 


4.1 










2/5 










2598.00 


4.6 


2601.25 


12.8 


/s + /e 














2671.52 


7.5 


/5 + /g 










2683.84 


17.3 


2685.94 


18.8 


/r, + /lo 










2693.91 


3.1 






/r. + /ii 










2697.29 


2.9 


2700.47 


5.6 


k + h 










2756.72 


3.7 


2756.21 


4.6 


h + /12 














2764.58 


3.8 


2/g 


2770.99 


1.7 


2771.25 


5.8 


2769.67 


2.3 


2770.64 


3.7 


/S + /l3 










2772.02 


3.0 






/g + /ii 










2783.12 


1.6 


2785.16 


3.0 


/6 + /l3 










2844.91 


2.0 






/g + /i3 


2859.62 


2.8 


2860.58 


2.1 










/5 + /l5 














2923.83 


4.7 


/4 + /16 










2929.43 


2.0 






h + fl6 










2932.71 


8.5 






fe + /16 










3005.60 


3.7 






/g + /16 










3008.70 


1.8 


3008.52 


4.8 


/g + /16 


3018.85 


12.3 


3019.23 


15.7 


3018.54 


10.4 






/ll + /l5 










3022.15 


3.3 


3023.05 


3.7 


/g + /i7 


3028.89 


3.5 


3028.58 


5.4 


3027.86 


2.0 






/l3 + fl6 


3107.48 


5.0 


3108.56 


4.1 












3117.52 


1.6 














/s + /18 










3186.34 


3.4 


3188.22 


5.0 


/l5 + /16 


3256.63 


3.1 


3257.12 


3.5 










2/16 


3266.71 


6.2 


3267.21 


6.2 


3267.42 


2.8 






/g + /18 


3272.97 


1.9 






3272.18 


3.2 


3272.91 


4.4 


/16 + fl7 


3276.75 


4.0 


3276.56 


3.7 










fl6 + /18 


3520.83 


4.2 


3521.39 


4.6 


3521.05 


4.2 






h + fl9 














3805.27 


2.8 


/g + /ig 


3890.36 


1.7 














/ll + /l9 










3903.26 


2.1 


3904.49 


3.1 


/ifi + /l9 


4138.22 


2.0 


4138.58 


2.8 










/l7 + /l9 


4148.26 


4.0 


4147.94 


3.4 










/18 + /l9 










4392.31 


1.0 







4.2 The second-order combination frequencies 

The combination frequencies appearing in Fourier analysis 
of the hght curves of pulsating white dwarf stars are most 
commonly interpreted as nonlinear distortions, most likely 
originating in their surface convection zones (Brickhill 1992). 



However, an alternative interpretation is resonantly coupled 
modes (Dziembowski 1982). These two hypotheses are of- 
ten difficult to separate, but as a general rule of thumb 
resonantly excited modes should have considerably larger 
amplitudes compared to signals that reflect light-curve dis- 
tortions. 
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Table 4. Third-order combination frequencies in our light curves of EC 14012— 1446 







200 1 


luiic 


200 1 


Ala\- 


2005 


V ) ■■ 

pii 


1 2007 


ID 


Frcq. 


Ampl. 


Frcq. 


Ampl. 


Frcq. 


Ampl. 


Frcq. 


Ampl. 




(mHz) 


(mmag) 


(mHz) 


(mmag) 


(mHz) 


(mmag) 


(A'Hz) 


(mmag) 


2/5 - /l7 














957.17 


3.4 


2/4 - /16 










957.74 


2.2 






h + k- fi6 










961.02 


4.3 






2/5 — /16 










964.30 


8.9 






2/b - /l5 










974.14 


3.8 






h + h — he 










1050.13 


3.3 






h + h — fi5 










1059.98 


2.5 






2/9 - fn 










1126.64 


2.7 






2/9 - /16 










1135.96 


3.5 






2/5 - /12 














1137.30 


4.8 


2/9 - /l5 










1145.81 


3.2 






/s + /e - /12 














1207.57 


4.4 


/4 + /5 - /9 










1209.89 


2.0 






2/5 - /9 










1213.17 


2.3 


1215.93 


4.9 


2/5 - fe 














1230.36 


8.9 


h + h — hi 










1285.55 


6.6 






h + /lo - /g 










1309.07 


6.1 






/s + /g - /e 














1315.05 


4.5 


2/16 - /18 


1379.24 


3.1 














2/g - h 










1470.67 


3.6 






/g + /16 - /i3 










1545.53 


2.1 






/s + /16 - /g 










1547.88 


2.7 






/) + /l7 - /l3 










1554.84 


2.3 






/ii + /ifi - /g 










1647.16 


6.0 






/g + /16 — /s 










1719.54 


1.7 






2/16 — /l3 


1792.59 


2.3 


1792.26 


2.4 










/i5 + /16 - /g 


1871.14 


2.1 














2/16 — /g 


1881.22 


5.0 


1881.58 


7.4 


1882.58 


1.8 






/16 + /i7 - /g 










1891.90 


3.2 






2/16 — /e 










1895.53 


2.0 






2/16 - /s 










1968.42 


0.8 






/16 + /ig — /g 










2753.84 


2.3 






/l + /l3 + /16 


3830.15 


3.2 


3830.15 


3.0 










/i + /g + /i7 


3850.16 


1.0 














3/5 














3901.87 


3.8 


2/5 + /g 










3982.84 


3.7 


3986.57 


7.7 


/s + 2/g 










4068.67 


3.9 


4071.26 


3.3 


2/5 + /13 














4074.29 


3.9 


/s + /g + /ii 










4082.12 


1.6 






/i + 2/16 


4087.98 


2.2 














2/5 + /16 










4231.71 


2.8 






/5 + /g + /16 










4317.55 


6.7 






h + ho + /16 










4327.62 


1.5 






2/9 + /16 


4404.35 


2.2 


4404.85 


5.3 


4403.38 


1.6 






/s + /l3 + /16 










4405.73 


1.5 






/g + /12 + /16 


4483.02 


0.8 














2/5 + /18 














4488.84 


2.5 


/"q -h f 1 -h /"l (S 

J y 1 J ±0 1 J -LU 


4492.98 


1.6 


4494.18 


2.3 










/s + 2/16 










4566.42 


1.6 






/s + /g + /18 










4571.18 


2.2 


4573.54 


3.7 


/g + 2/16 


4652.21 


4.9 


4652.83 


5.1 


4652.25 


2.0 






/g + /16 + /i7 


4662.25 


1.5 


4662.19 


2.1 










/i + /g + /ig 


4711.62 


0.9 














/is + 2/16 


4740.83 


1.7 














/s + /16 + /18 










4820.05 


1.5 






3/16 


4900.07 


1.9 














/g + /16 + /18 










4905.88 


2.0 






/l3 + /16 + /l7 


5473.54 


2.6 














h + hi; + /i<) 






5524.21 


1.8 










/ifi + /i7 + /i:j 


oV81.(il 
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Table 5. Fourth and fifth-order combination frequencies in our light curves of EC 14012—1446 



April 2001 June 2001 May 200r) April 2007 



ID 


Frcq. 


Ainpl. Frcq. 


Ampl. 


Frcq. 


Ampl. 


Frcq. 


Ampl. 




(mHz) 


(mmag) (AtHz) 


(mmag) 


(mHz) 


(mmag) 


(mHz) 


(mmag) 


/n — P Fir; — 2. T A 1 




69.32 


5.5 










2 ft, — f r — A 1 








72.38 


1.8 






Tfi — P T^ £i — ?. IK 








420.54 


2.6 














669.41 


1.8 






2fifi — f A — ft; 

iO J 4 JO 








672.69 


2.1 








2216.80 


2.1 












3 ft; — fifi 








2263.30 


2.7 






f. _p f_ _p f„ _ f,,, 

J4 Jo J9 J 10 








2345 85 


2.2 






/. TK — P TCi — /I f: 

•^joT^jy jro 








2349.13 


3.3 






ft; + fo + fu — flfi 

Jo^Jo^jy Jib 








2432.10 


2.1 






fr; + 2 fa — f 1 fi 
J i> ^ J y J ID 








2434.97 


2.0 






2 fc; -P fn — ffi 












2615.68 


3.2 














2853.56 


3.1 


y J lo JO 








2943.69 


2.2 






/4 + 2/i6 — /g 








3178.31 


1.2 






3/l6 — h 


3514.57 


2.7 3515.19 


4.3 










fl + 2/9 + /l6 


5225.61 


1.7 












3/5 + /9 












5287.19 


2.7 


2/5 + 2/9 








5367.68 


1.8 






2/5 + /9 + /16 








5616.55 


2.5 






/s + 2/9 + /16 








5702.38 


1.9 






3/9 + fid 


5789.85 


1.2 












fr, + fn + 2flfi 








5951.26 


1.0 






/r, + 2/9 + /18 








5956.01 


1.0 






2/9 + 2/16 




6038.46 


3.0 










h + h + fi6 + /18 








6204.89 


1.0 






3/16 + /9 


6285.57 


2.0 












/9 + 2/13 + /16 


6689.76 


1.1 












fg + fi3 + /16 + fi7 


6859.04 


1.4 












2/5 + /9 - /ll - /16 








950.85 


2.6 






3/9 — /s — /18 








968.17 


1.8 






3/5 + /9 - /16 








3648.13 


2.3 






/s + /ll + /16 + /l9 — /18 








4948.63 


1.3 






3/5 + 2/9 








6666.68 


1.1 






3/5 + /9 + /16 








6915.55 


0.8 






2/6 + 2/9 + /16 








7001.39 


1.1 







Figure 4 shows the ratio of the amphtudcs of the indi- 
vidual combination frequencies with respect to the prod- 
uct of the amplitudes of the parent signals. Note that 
these amplitude ratios are given in units of mag~^, whereas 
the amplitudes themselves are quoted in millimagnitudes. 
Whereas 85% of the data points axe located in an interval 
of 8.5 ± 6.2 mag~^, the remaining combinations show con- 
siderably larger amplitude ratios. 

By far the largest amplitude ratio is duo to the com- 
bination frequency at 2196 /iHz present in both data sets 
from the year 2004. It therefore is our best candidate for a 
resonantly coupled mode. In addition, the second-order com- 
bination at 2296 /iHz also stands out in Fig. 4. Both of these 
combinations include the lowest-frequency independent sig- 
nal fl = 821 /iHz. This latter signal also has two other com- 
binations with rather large amplitude ratio with respect to 
its parent variations (at 2207 and 2455 /iHz. respectively). 
The third combination that shows enhanced amplitude in 
both data sets from 2004 is the one at highest frequency, 
involving /19. Finally, the combination frequency diflference 



at 268 /iHz in June 2004 also has large amplitude, but the 
residual amplitude spectrum in this data set has high noise 
around this frequency. We therefore suspect that the large 
amplitude of the 268 /iHz signal may be spurious. 

Turning to the data sets from 2005 and 2007, we do 
not find such extreme amplitude ratios. Three of the five 
high-amplitude combination frequencies again involve /19, 
and the other two contain "satellite" frequencies around the 
strong signal at /9 « 1385 /iHz. Of highest interest is the 
combination at 1133.75 /tHz in the April 2007 data. Its fre- 
quency is consistent within the errors with that of the in- 
dependent signal /2 that was only present in the year 2004. 
However, in 2004 only one of the parents of the combination 
present in 2007 was detected. Therefore f'z must be regarded 
as an independent signal in the data from 2004. 

We remind that during the four epochs of observations 
of EC 14012—1446 at our disposal, amplitude variations oc- 
curred. This puts us into a position to examine whether 
or not the amplitudes of combination signals present in 
more than one of these data sets followed the corresponding 
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Figure 3. Amplitude and frequency variations of six independent signals in the light curve of EC 14012—1446 during the observations 
in May 2005. The error bars of the individual data points are la. The panels are grouped with respect to the signal, and have been 
arranged to compare the amplitude and phase behaviour directly. Note the different ordinate scales in the different panels. The central 
subset of data used for frequency analysis corresponds to 4 < Time < 10 in this graph. 



changes of their parent signals. We restrict this analysis to 
the combination frequency sums because these are less af- 
fected by observational noise, i.e. have better signal-to-noise 

on average. 

In order not to make any assumptions about the origin 
of the combination frequencies, we chose to compare the 
individual amplitude ratios of the combinations that occur 
in more than one data set with respect to their mean and 
their deviations from their formal accuracies. If there were 
no changes in the relative amplitudes, we would expect a 
normal distribution to result. The outcome of this test can 
be found in Fig. 5. 

Most of the amplitude ratios follow a Gaussian distri- 
bution with a width of fs 1.2a. Given the fact that formal 
errors of the determinations of amplitudes tend to underes- 
timate the real errors, this is not surprising. 



There are nine recurring second-order combination sig- 
nals whose relative amplitude varies significantly, and they 
correspond to four different sets of parents: /s -|- /g, /e + /g, 
2/g , and /g -I- fie . All these combinations contain the mode 
/g, but we failed to find any systomatics in the temporal 
behaviour of the relative amplitudes of its combinations due 
to the sparsity of data. 

We also tested whether possible interactions with other 
combinations that happen to have frequencies close to the 
four signals under investigation might be able to affect their 
amplitudes. The result of this test was negative because in 
general all possible alternative parent frequency combina- 
tions have much lower amplitudes than the parents we as- 
signed. The single exception to this rule is 2/g in the 2007 
data set, whose amplitude could be mildly affected by inter- 
ference with fe + fii- 
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Figure 4. Amplitude ratios between the second-order combination frequencies and the product of their parent signal amplitudes, with 
1(7 error bars indicated. Multiple determinations at a single frequency are due to combinations present in more than one data set. 
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Figure 5. The distribution of the deviation of the amplitude 
ratios of recurring second-order combination frequencies with re- 
spect to their mean and formal errors (histogram bars). For com- 
parison, two Gaussians with standard deviations of Icr (dotted 
line) and 1.5(t (dash-dotted line) are shown. 



4.3 The third-order combination frequencies 

The third-ordor combination frequencies determined can be 
separated into two approximately equally populated groups: 
combinations of the form /a + /s + fc and fA + fB — fc, 
respectively. An important fact about the second type of 
combination signals it that their frequencies are in the same 
range as those of the independent signals (see also Vuille et 
al. 2000). Consequently, they may be mistaken for indepen- 
dent frequencies if not systematically searched for. 

As a matter of fact, such confusion has occurred in 
the preliminary report on EC 14012—1446 by Handler & 
Romero- Colmenero (2005), where two signals now identified 
as third-order combination frequencies were believed to be, 
and interpreted as, independent mode frequencies. We illus- 
trate the severity of this problem in Fig. 6, that shows the 
schematic amplitude spectrum of our multisite data set with 
the independent frequencies and combination signals sepa- 
rated. Most of the signals in this frequency range in fact 
correspond to combination frequencies! 



We again stress that we prefer to misinterpret a peaJi in 
the amplitude spectrum as a combination signal as opposed 
to falsely taking combination frequencies as real modes; we 
followed this guideline throughout the frequency analysis. 
We will discuss the astrophysical implications of confusing 
combination signals with independent mode frequencies be- 
low. 



4.4 The fourth and fifth-order combination 
frequencies 

We were able to identify several high-order combination fre- 
quencies in the data. To our knowledge, this is the first re- 
port of the detection of fifth-order combination signals in the 
light curves of a pulsating white dwarf star. The identifica- 
tion of the signals that are pure frequency sums should be 
secure because there are no ambiguities in assigning the par- 
ents, whereas some identifications involving frequency sums 
and difl^crences may be incorrect. However, it is important 
to point out that we only matched significant peaks in the 
amplitude spectrum, and that we succeeded in explaining 
them all. 



5 INTERPRETATION 

5.1 The independent frequencies 

With our approach to frequency analysis, most of the inde- 
pendent frequencies should correspond to normal mode os- 
cillations. We interpret /i, /2, /is and /19 as single modes. 
The triplets /i2,/i3,/i4 and /i5,/i6,/i7 are normal mode 
multiplots; their splitting of « 10 ^Hz is most likely duo to 
rotation. We believe that frequency /g is also caused a nor- 
mal mode. It does show a ~ 10 /iHz splitting at times, but 
from time to time a « 14.5 fiHz splitting is present. The sig- 
nal /s is also suggested to be caused by a normal mode, and 
it once showed a companion frequency again separated by 
« lO/jEz. 

Assuming that the modes where we found rotational 
splitting have a spherical degree £ = 1, we can estimate the 
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Figure 6. Schematic amplitude spectrum of the May 2005 data (see Tables 2 and 4). The independent frequencies are the thick full 
lines, whereas the third-order combination frequencies are shown with thin dotted lines. 



rotation rate of the; star. The average amount of the split- 
ting we believe to be of rotational origin is 9.94 ± 0.12 pHz. 
Taking into account that this is about 4% off the asymp- 
totic value (Brassard ot al. 1992) and that the rotational 
splitting is somewhat variable, we arrive at a rotation fre- 
quency of 19.1 ± 0.8/LtHz or at a fairly short rotation period 
of 0.61 ± 0.03 d. 

Some of the independent signals in our light curves 
(/4) /fi, fs, fii) did not conform to rotationally split peaks as 
interpreted before. There may be two possibilities for their 
origin. Firstly, they could be members of rotationally split 
multiplets with different £. However, no reasonable assump- 
tion (either the ~ 10 /iHz splitting is due to £ = 1 modes 
or it is due to the i = 2 rotational splitting) in this frame- 
work can explain the whole set of signals detected around /g . 
Secondly, these unexplained signals may be artefacts due to 
amplitude and/or frequency variations. Wo recall from Sect. 
4.1.1 that although such temporal variability in the ampli- 
tude spectra is definitely present, we cannot prove this hy- 
pothesis. In addition, the presence of a "pure" triplet with 
a 14.5 /iHz spacing in April 2007 is hard to explain with 
amplitude/frequency variability only. 

The period spacing between the normal modes of a pul- 
sating white dwarf star serves to obtain its mass if it allows a 
determination of the asymptotic g-mode period separation. 
Following the discussion above, we interpret /i, /2, /s, /g, 
/i3i /i6, /i8 and /ig as normal mode frequencies with an 
azimuthal order m — 0, but with different radial overtone 
k and/or spherical degree I. The assignment of m follows 
from the multiplct structure around /g, /13 and /ig (in case 
they are £ = 1 modes), but is assumed for the others. How- 
ever, if wrong, this assumption would not significantly affect 
any mass determination, because the frequency difference of 
consecutive radial overtones is significantly larger than the 
rotational splitting. 

In any case, our attempts to establish the asymptotic 
g-modc period separation of EC 14012—1446 failed. There 
may be two possible causes: firstly, the effects of mode trap- 
ping (e.g. see Brassard et al. 1992) could be so strong that no 
regular period spacing can be found within the few modes 
available. Secondly, not all of the modes we adopted may 



Table 6. Changes in the light curves of EC 14012-1446. 



Month/ Year 


WMP 




rms scatter 




(s) 


(mmag) 


(mmag) 


June 1994 


598.4 


66.9 




April 2004 


658.9 


59.9 


48.9 


June 2004 


663.6 


64.9 


53.6 


May 2005 


684.4 


72.2 


60.3 


April 2007 


692.0 


85.7 


71.2 



have the same spherical degree i. A comparison of the pe- 
riods of these eight modes with the ones identified as f = 1 
by Kleinman et al. (1998) in G 29-38 shows good agree- 
ment except for /ig, assuming the stars arc homologous. 
We conclude that an asteroseismic mass determination for 
EC 14012—1446 must either await the detection of a more 
complete set of normal modes or requires systematic match- 
ing with a model grid (Metcalfe, Montgomery & Kanaan 
2004, Castanheira & Kepler 2008). 

With the independent frequencies and their amplitudes 
determined, we can compute the weighted mean pulsation 
period of EC 14012—1446 at the different observing epochs, 
adopting the definition by Mukadam et al. (2006): 

WMP= (1) 

where WMP is the weighted mean period, the Pi are the pe- 
riods of the individual pulsation modes and the Ai are their 
respective amplitudes. We also computed the rms scatter 
of each individual data point in our light curves and the 
square root of the total power (X)(^0^)^^^ ^ indications 
of the pulsation power of the star. The first measure makes 
the (justified) assumption that the intrinsic variations in our 
light curves dominate the measurement noise and it refiects 
the total variability. The second measure requires knowledge 
of the independent mode frequencies, but allows a direct 
comparison with the work by Mukadam et al. (2006). 

The results of these calculations are given in Table 6. 
We also computed the WMP and the square root of the total 
power from the data by Stobie et al. (1995) for complete- 
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ness. As these authors quoted nightly amphtudes for their 
independent frequencies, wo adopted the mean of these am- 
phtudes to obtain the WMP and {J^i'^if V^^ this ob- 
serving period. 

Tabic 6 shows a clear trend: the mean period of the 
variability increased from one observing epoch to the next 
(a fact that can already be seen by eye in Fig. 2). The to- 
tal light amplitude increased during our observations, but 
was at an intermediate level when Stobie et al. (1995) ob- 
served the star. As with the apparent frequency changes of 
some individual modes (Sect. 4.1) we caution against inter- 
preting the change in the WMP as a monotonous trend in 
time: the gaps between our observing runs were considerably 
larger than the observing runs themselves. Because we have 
already found amplitude and frequency changes within one 
observing run, it is quite possible that the star altered its 
behaviour on time scales not sampled by us. 

If we compare these results with the sample of white 
dwarfs analysed by Mukadam et al. (2006), we find that 
EC 14012-1446 is among the highest-amplitude ZZ Ceti 
stars and that it is located close to the centre of the insta- 
bility strip, consistent with spectroscopic results (Fontaine 
et al. 2003, Bergeron et al. 2004). The high amplitudes let us 
speculate that most pulsation modes are indeed £ = 1, suf- 
fering the smallest geometrical cancellation. The temporal 
variation of the WMP we observed is not unusual compared 
to other ZZ Ceti stars (Mukadam et al. 2007). 



5.2 The second-order combination frequencies 

Most of the second-order combination signals have ampli- 
tudes that scale with the product of the amplitudes of their 
parent modes. These are consistent with an interpretation in 
terms of light-curve distortions. Consequently, we can try to 
determine some convection parameters with light-curve fit- 
ting following Handler, Romero-Colmenero & Montgomery 
(2002), based on Wu's (2001) analytical expressions for the 
combination signals' amplitude and phases (also see Yeates 
et al. 2005). 

Using the data from April 2004, assuming all parent 
modes are i = 1 and only using the mode triplets, we 
find a response parameter 2/3 + ^ = 9.8 ± 0.3, convective 
thermal time Tc = 135 ± 10 s, and an inclination angle 
e = 29.5 ± 1.5° (see Wu 2001 for the definitions of these 
parameters) when fitting the light curve. For 2005 however, 
no satisfactory fit could be obtained. 

As noted before, our best candidate for a resonantly 
coupled mode is the 2196 /iHz signal in the 2004 data, as 
its relative amplitude is by far the highest, and since it var- 
ied in accordance with the amplitudes of its parents /i and 
fr. However, the interpretation of this 2196 /xHz signal may 
be a chicken-and-egg problem: in June 2004, this signal had 
about 70 per cent of the amplitude of its parent mode /i. 
Given that this number is close to unity and that different 
spherical degrees, and thus different geometrical cancellar 
tion, of the parent and the combination frequencies, may 
come into play, some doubt as to which signal is indeed the 
parent and which one is the combination remains. In addi- 
tion, a 2196 /xHz signal fits much better into the sequence 
of the observed mode frequencies of EC 14012—1446 than 
/i = 821 /iHz does. 



5.3 The third-order combination frequencies 

The large rmmber of third-order combination frequencies de- 
picted in Fig. 6 follows a simple structure: all these groups 
of peaks are separated from the dominant modes by the fre- 
quency differences between the three strongest modes and 
their rotationally split components. This means we see a 
pattern fairly regularly spaced in frequency, which, to first 
order, is also a pattern regularly spaced in period. 

If those third-order combinations are not correctly iden- 
tified, but are misinterpreted eis independent frequencies, an 
incorrect determination of the period spacing of the nor- 
mal modes may therefore result. As an experiment, we de- 
termined the mean period spacing of all independent and 
third-order combination frequencies detected between 950 
and 2550 /itHz in the May 2005 data. We arrived at a for- 
mally significant mean period spacing of 17.7 s. Some "radial 
overtones" in this scenario would be missing, but no inte- 
ger multiple of 17.7 s would give a reasonable fit. Inciden- 
tally, this completely incorrect period spacing is close to the 
true period spacing of the massive ZZ Ceti star BPM 37093 
(Kanaan et al. 2005). Hence, EC 14012-1446 could be mis- 
interpreted as having high mass, inconsistent with spectro- 
scopic results (Fontaine et al. 2003). It is therefore essential 
that combination frequencies are correctly recognized in the 
amplitude spectra of pulsating white dwarf stars. 

Vuille et al. (2000) pointed out a similar problem in 
their study of the prototype pulsating DB white dwarf star 
CD 358. However, there are some difi^erences between our 
results and those by Vuille et al. (2000): the third-order 
combination frequencies revealed in their work only man- 
ifested themselves as "odd" peaks lying close to or within 
rotationally split multiplets, and these authors could explain 
all of these "odd" peaks with combination frequencies. In 
our case, the third-order combinations create apparent ad- 
ditional frequency multiplets, and we found signals within 
the rotationally split normal mode multiplets that we were 
unable to explain by combination frequencies; we suspect 
they axe artefacts caused by amplitude/frequency variations 
during the measurements. 



6 CONCLUSIONS 

We have acquired more than 200 hours of time-resolved pho- 
tometry of the ZZ Ceti star EC 14012-1446. The star turned 
out to have a complex pulsation spectrum, that we were able 
to understand to some extent, but not to our full satisfac- 
tion. 

We detected eight normal modes of the star with differ- 
ent radial overtones. Not all of them must necessarily have 
the same spherical degree £. Some of the normal modes are 
arranged in groups, in most cases showing separations of 
~ 10 /iHz. We interpreted that as the sign of rotational m- 
mode splitting, and inferred the rotation period of the star 
(0.61 ± 0.03 d), assuming modes oi £ = 1. 

We were unable to determine a regular period spacing 
between the pulsation modes of EC 14012—1446, and there- 
fore cannot provide a mass estimate. The reason may be 
mode trapping or that we see a mixture of modes with dif- 
ferent spherical degree £. 

The pulsation spectrum of the star changes over time. 
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We have measured changes of ±8% in the weighted mean 
period and of ±20% in the square root of the total pulsa- 
tion power. Our data have insufficient sampling and time 
base to determine the time scale of the variations in the 
normal modes. However, wc found systcmatics in the tem- 
poral changes of the pulsation spectrum: the higher the ra- 
dial overtone of a mode, the more pronounced the variations 
are and amplitude and frequency changes go hand in hand. 
The dependence on radial overtone is easy to understand: 
since the kinetic energy of the g modes in a pulsating white 
dwarf star becomes weighted towards the outer regions with 
increasing radial overtone, less energy is required to disturb 
a high-overtone mode. 

Most of the second-order combination frequencies have 
amplitudes consistent with the hypothesis that they origi- 
nate from light curve distortions of the star. However, some 
could be resonantly excited modes because they have rel- 
ative amplitudes much higher than the other second-order 
combinations. We even found one case where it was no longer 
clear which was the parent mode and which was the combi- 
nation, although this could partly be an effect of geometrical 
cancellation. 

The amplitudes of the second-order combinations varied 
in time in accordance with their parent modes. The excep- 
tions to this rule are a number of combinations that involve 
the mode /g. This mode is also special in the sense that it 
shows complicated "multiplet" structure in one of our data 
sets, that wc tentatively attributed to the effects of ampli- 
tude and frequency variations during the observations. 

Finally, we showed how a misinterpretation of third- 
order combination frequencies of the form fA + fs — fc as 
normal modes, can lead to a completely nonsensical deter- 
mination of the stellar mass. Combination frequencies in- 
truding into the domain of the normal modes in the Fourier 
spectra of the light curves must be carefully searched for, 
and eliminated from an asteroseismic analysis. 

The ZZ Ceti star EC 14012-1446 remains an attractive 
target for further study, as we could only partly understand 
its pulsations. The high pulsation amplitudes of the star may 
facilitate the application of Montgomery's (2005) method 
to determine convection parameters and mode identifica- 
tions from high signal-to-noise light curves. The detection 
of more normal modes of pulsation and the correct identi- 
fications of at least some of them may permit a full astero- 
seismic analysis. A reliable determination of the time scales 
of the amplitude and frequency variations of the star would 
be worthwhile to pin down their origin, but would require a 
data set with a time base in excess of two weeks, excellent 
duty cycle at the same time, as well as a co-operating tar- 
get. Under such circumstances, precision asteroseismology 
of a ZZ Ceti star may be feasible. 
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